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ABSTRACT 

Several scenarios have been suggested to explain the phase-space distribution of the 
Milky Way (MW) satellite galaxies in a disc of satellites (DoS). To quantitatively 
compare these different possibilities, a new method analysing angular momentum di- 
rections in modelled data is presented. It determines how likely it is to find sets of 
angular momenta as concentrated and as close to a polar orientation as is observed for 
the MW satellite orbital poles. The method can be easily applied to orbital pole data 
from different models. The observed distribution of satellite orbital poles is compared 
to published angular momentum directions of subhalos derived from six cosmological 
state-of-the-art simulations in the Aquarius project. This tests the possibility that 
filamentary accretion might be able to naturally explain the satellite orbits within the 
DoS. For the most likely alignment of main halo and MW disc spin, the probability 
to reproduce the MW satellite orbital pole properties turns out to be less than 0.5 
per cent in Aquarius models. Even an isotropic distribution of angular momenta has a 
higher likelihood to produce the observed distribution. The two Via Lactea cosmolog- 
ical simulations give results similar to the Aquarius simulations. Comparing instead 
with numerical models of galaxy-interactions gives a probability of up to 90 per cent 
for some models to draw the observed distribution of orbital poles from the angular 
momenta of tidal debris. This indicates that the formation as tidal dwarf galaxies 
in a single encounter is a viable, if not the only, process to explain the phase-space 
distribution of the MW satellite galaxies. 

Key words: galaxies: kinematics and dynamics - galaxies: dwarf- galaxies: formation 
- Local Group - dark matter 



1 INTRODUCTION 

The satellite galaxies of the Milky Way (MW) are dis- 
tributed in a highly inclined plane a round the MW disc 
l|Lvnden-Belllll976l : iKroupa et aLlbOOSl ). termed the Disc of 
Satellites (DoS) . This is tru e individually for the 11 'clas- 
sical' satellites l|Metz et al.l i2007) as well as for the fainter 
ones detected in the Sloan Digital Sky Survey (|Kroupa et al.l 
I2OIOI ). In addition, globular clusters of the MW categorized 
as young halo clusters populate the same plane, and streams 
of stars and gas show a preference to align with it, too, thus 
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being evi dence for a vast polar structure (VPOS) around 
the MW (jPawlowski et al.ll2012l ). 

This strong spatial anisotropy is supported by the mo- 
tions of th ese satellite galaxi es. Using proper motion mea- 
surements, iMetz et al.l ( 20081 . M08 hereafter) have derived 
the orbital poles (direction of angular momenta) for eight 
satellite galaxies. They found a strong alignment of six or- 
bital poles close to the normal vector to the DoS, indicating 
that it is rotationally supported. In addition, the Sculptor 
dwarf galaxy is counter-orbiting, but also within the DoS. 
Thus, seven of eight satellites with measured proper motions 
have aligned orbital axes. 

Several attempts have been made to explain the afore- 
mentioned spatial and orbital anisotropy within the stan- 
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dard model of cosmology. Dwarf Ralaxies miRht be acc reted 



in groups l|Li fc HelmfeoOSi : iD'Onghia fc Lakell2009l ), but 
iMetz et al.l (l2009l l have shown that such groups, observed 
extra-Galactic dwarf galaxy associations, are far too ex- 
tended to explain the thin DoS structure. Comparing the 
orbital energy of MW satellites with dark matter su bhalos 
from the Via Lactea II simulation. iRocha et al.l ('2011) find a 
wide spread in satellite infall times. This is in conflict with a 
collective accretion of the satellites in a single group. A com- 
mon infall of MW satelli tes as former satellit e galaxies of the 
LMC is investigated bv lNichols et all l|201ll ). who conclude 
that "the extended disk-of-satellites cannot be explained by 
the dwarfs being bound to the LMC within the last two 
apogalacticons, and may have another origin." Therefore, 
the idea of a group of galaxies centred on the LMC which 
was accreted is strongly disfavoured. 

iLibeskind eFal] (|2009l ') demonstrate that it is in prin- 
ciple possible that dark matter (DM) dominated satellites 
are similarly aligned as the most luminous MW satellites. 
Howe ver, this was s hown to be an unsatisfactory solution, 
since iKroupa et all (|2010l ') demonstrated that their results 
indicate that only 0.4 per cent of all existing cold dark mat- 
ter (CDM) halos of MW mass would host a galaxy simi- 
lar to the MW wi th a similar spat i al dis tribution of satel- 
lites. Furthermore. iLlbeskind et al.l (|2009l ) only resolve very 
massive subhalos (over 2.6 x 10^ Mq), which is inconsistent 
with the lower dynamical mass estima t es of most MW sate l- 
lites (|Mateolll998l : IWalker et al]|2007l : [Strigari et al.ll20oi ). 
The significant over-abundance of predicted bright satel- 
lites compared t o obse rvations has first been identified by 
iBovill fc Ricottil (|201ll ') and termed the "bright sateUite 
probl em". In addition to this discrepancy in absolute num- 
bers, iKroupa et ahl l|2010l ) have shown the inconsistency of 
the predicted mass-functions of luminous subhalos with that 
of observed MW satellite galaxies. 

A series of logically incompatible conclusions have 
been reached by authors working in the ACDM frame- 
work att empting to explain the M W satellite population: 
iNichols fc Bland-Hawthorq (120111 ) argue that the dSph 
satellites had to have fallen in at a redshift z > 3 in order for 
their gas content to be r emoved by means of tidal and ram- 
pressure stripping, while [Peason et al] (|201ll ') argue that the 
same dSph satellites must have fallen in recently in order 
for them to form the DoS and they discuss an example of 
group infall at z — 0.6. Apart from these mutually exclusive 
results, it has transpir ed that the proper motions of dSph 
satellites exclude infall (|Angus et al.ll201ll ) . From Via Lactea 
II data Hensler fc Petrov (2012, in preparation) model the 
subhalo system around the MW and demonstrate that all 
DM subhalos should experience star formation, do not suffer 
tidal disruption, but only ram-pressure gas loss, and survive 
as faint dSphs. The most important issues are that all satel- 
lite dSph galaxies with subhalo masses larger than 10^ Mq 
have present-day surface brightnesses detectable by SDSS 
data and are isotropically distributed. 

In one of the most recent attemp ts at explaining the 
DoS within ACDM. iLovell et all (|201ll . Lll hereafter) have 
reported that CDM simulations naturally lead to satellites 
coherently rotating in quasi-planar distributions, such as 
the satellites of the MW. They calculated the directions 
of orbital angular momenta of DM subhalos in six galac- 
tic high-resolution halos taken from the Aquarius project 



(|Springel et ahllioOSl ). These are CDM simulations of similar 
halo masses to the assumed DM halo of the MW. Six differ- 
ent models have been analysed at resolution level 2, called 
Aq-A2 to Aq-F2. For details on the cosmological simulations 
and the method of determining subhalo angular momenta, 
the reader is referred to Lll. In that paper, while using 
the results to interpret the MW data, a quantitative test 
showing whether or not the angular momenta of the subha- 
los could feasibly represent the distribution of MW satellites 
was not given. Here we make that test, in particular because 
the Lll results have been used to ar gue that the DoS is nat- 
urally explained within ACDM (e.g. lKeller et al.ll2012l '). 

Motivated by the DoS and the apparently coherent 
orbits of the satellite galaxies within it, a n alternative 
scena rio of their origin has been pro posed ('Lvnden-Belll 
1 19761 : iKroupa et al] l2005l : iMetz et al.1 2008; Kroupa et al] 
[2OI0I). They might be tidal dwarf galaxies (TDGs), formed 
from the tidal debris in an encounter between t he early 
MW and another, still gas-rich (proto-) galaxy (|Zwickvl 
'1956'). Galaxy collisions can lead to perpendicularly ori- 
ented discs of debris, with aligned o rbits as well as counter- 
orbiting material as demonstrated in lPawlowski et ah l (|201lL 
Pll hereafter). Such encounters are observed even in the 
prese nt epoch, for exam ple in the interacting system VV 
340 (lArmus et al.l bOOgl). TDGs are found to form both 



in simulated dWetzstein et al.ll2007l : iBournaud fc Duj 



Bournaud et al.l l2008t) and observed (iMirabel et aT 



2006 



1992 



Hunsberger et al.l 19961 : IWeilbacher et al.l 2003f) galaxy col- 
lisions. They seern to be long-lived obje c ts jlCroupal 1997 



iRecchi et al] l2007l : ICalianni et"al] I2OI0I : IDuc et al.l I2OII ) 

and share the same properties as dwarf elliptical galaxies 
(Dabringhausen et al., in prep.). This makes them a fun- 
damental addition, if not even an alternative, to cosmo- 
logically formed dwarf galaxies. In fact, that TDGs may 
be the dominant satellite dwarf galaxy population should 
be a necessary outcome o f ACDM structure formation 
(IQkazaki fc Taniguchillioool ). 

Here, a new method is presented which determines how 
likely it is to find as strongly a clustered distribution of or- 
bital angular momentum vectors similarly close to the equa- 
tor of the MW as the MW satellite orbital poles, if the angu- 
lar momenta of the satellite population were drawn from the 
distributions in given models. This method is then applied 
to the angular momentum directions of ACDM subhalos as 
determined and published by Lll, to those derived from 
the two Via Lactea simulations (|Diemand et al. I l2007l . l2008l ) 
and to those of tidal debris as determined from models of 
galaxy interactions. Section[2]describes the assumptions, the 
observed situation, the models and the method. In Section 
[3] the results are presented, followed by the conclusions in 
Section ID 



2 ANALYSIS 

2.1 Assumed galactic disc orientation in a DM 
halo 

It is generally assumed that the angular momentum of the 
baryonic disc galaxy aligns with the angular momentum of 
the DM halo in which it resides. This assumption seems to 
be supported by the alignment of the minor axis of the inner 
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Figure 1. The data of figure 1 in Lll, the cosine of the angle 
between the main halo spin and the subhalo angular momentum 
vector ©H Si are plotted here as a histogram (input data courtesy 
of Mark R. Lovell). For this plot, the subhalo orbital angular mo- 
menta of the six individual Aquarius simulations of Lll have been 
co-added. This results in the black, solid histogram, which gives 
the average probability density of the six simulations. The varia- 
tion of the different simulations is illustrated by the grey, shaded 
area. For each bin, it illustrates the maximum and the minimum 
of the probability density determined for the six Aquarius halos 
individually. Overall, the probability density of the simulations 
lies close to the isotropic case, illustrated by the horizontal dotted 
line. Co-rotating orbits are slightly preferred in all simulations. 
The deviation from the isotropic value is strongest for the bins 
closest to the DM spin axis, where the difference in probability 
density is up to a factor of two. Thus, while there is an over- 
abundance of subhalos with spins close to the main halo spin, all 
other spin directions are present in significant fractions. Note that 
the sum over all angles (of the area below the histogram) yields 
a probability of 1. 



part of DM halos with the axis of the dis c galaxies in them 
ijBailin et al.ll2005h . lSharma fc Steinmetj f2005l 1 find a good 
correlation between baryonic gas and DM halo spin, with a 
mean misalingmcnt-angle of only 18.9° at a redshift of 0. The 
work of03ctt ct al. (2010) also shows that in simulations it is 
most likely for a galaxy to have its spin vector aligned with 
the direction of angular momentum of its parent DM halo. 
In particular the angular momentum of the inner halo (0.25 
times the virial radius) aligns well with the galaxy spin, the 
median angle of misalignment being 23.9°. Taking the angu- 
lar momentum of the whole halo, the median misalignment 
rises to 34.4° , but still the majority of galaxy spins are close 
to the main halo spins. Perpendicular orientations are rare 
(12 per cent within ±15° arou nd the perp endicular orienta- 
tion according to figure 17 in lBett et al.l|2010 ), just as are 
dark halos in which the inner (^ 0.25 rvir) and outer regions 
are tihed by ^ 90°. 

In the following, it is therefore assumed that the galaxy 
rotates in the same direction as the parent halo, in line with 
what Lll state in their paper. It is assumed that the MW 



disc spin and the main halo spin in the simulations are par- 
allel. In Section r3.4l the requirement of this most likely align- 
ment will be dropped, instead choosing an orientation of the 
disc galaxy which is most favourable to the formation of po- 
lar orbits, but puts the main halos spin and disc galaxy spin 
at an unlikely 90°-angle. 



2.2 The MW satellite orbital poles 

The distribution of the MW satellite orbital poles (M08) 
shows two characteristic properties. First of all, the orbital 
poles cluster close to the normal to the MW DoS. The six 
best-aligned orbital poles show a spherical standard distance 
(M08) of only A^JJ' = 35.4°. Secondly, the orbital poles 
preferentially fall close to the MW equator, the satellites 
thus move along on polar orbits. This is obvious from the 
direction of the average orbital pole for the six best-aligned 
orbital poles, it has an angular distance of only dMW ~ 9.4° 
from the MW equator. 

New proper-motion measurements for the satellite 
galaxies might lead to updated orbital pole directions, which 
in turn might change these values. The uncertainties of the 
orbital poles are the projected uncertainties of the angu- 
lar momenta of the satellite galaxies. They are dominated 
by the, often large, uncertainties in the measured proper 
motions. M08 have shown that the uncertainties in the dis- 
tance, position and radial velocity of satellite galaxies, the 
distance of the Sun from the Galactic centre and the circu- 
lar velocity of the local standard of rest are negligible com- 
pared to the infiuence of the proper motion uncertainties on 
the directions of the orbital poles. But large proper-motion 
uncertainties do not necessarily result in large orbital-pole 
uncertainties. What is important for the determination of 
the orbital pole of a satellite galaxy is the direction of its 
motion with respect to its position relative to the Galactic 
centre, not its velocity. 

As a result, the orbital poles of most satellites co- 
orbiting in the DoS are relatively well defined (directional 
uncertainties < 15°). Only Draco and Carina show large un- 
certainties. However, their uncertainties fall on great-circles 
passing close to the average orbital pole of these six best- 
aligned orbital poles (see figure 1 of M08). This is why, on 
average, varying their orbital poles within the uncertainties 
cancels out when determining and diiiwH- Within the 

current orbital pole uncertainties, updated proper motions 
can result in both larger or smaller values for A^^ and 
dMW, but the changes can be expected to be small. There- 
fore, in the following analysis the given values for the two 
parameters are adopted. 

Figure [1] is a remake of figure 1 in Lll using their 
data. It plots the distribution of the cosine of the angle be- 
tween the main halo spin and the orbital angular momenta 
of individual subhalos, cosQh s- All models peak close to 
cosQh-3 = 1, so it is most likely for a subhalo to orbit in 
the same direction as the main halo. Three of the models 

^ Note that this is not inevitable. If the uncertainties would have 
been oriented perpendicular to this orientation, a variation of the 
satellite orbital poles along the uncertainties would preferentially 
lead to larger A^^-values as the distance between the orbital 
poles would preferentially increase. 
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Figure 2. The same plot as shown in Figure [T] but note the 
different range on the Y-axis, fn addition, this plot includes the 
orbital poles of the eight MW satellite galaxies (dashed black 
histograms with error bars) from M08, generated in the same 
way as the model histograms, fn the upper panel it is assumed 
that the spin direction of the Milky Way galaxy, which points 
to the southern Galactic pole, aligns with the main halo spin, 
as Lf f have argued. Therefore, the histogram shows the distance 
of the MW satellite orbital poles measured from the MW south 
pole. The distribution found in the simulations clearly differs from 
the observed one of the MW satellites. The MW satellites are 
preferentially on polar orbits, perpendicular to the MW disc spin 
and by implication to the spin of its supposed DM halo, at odds 
with the Lll results. Most of the eight known satellite orbital 
poles point close to the MW equator, about 90° away from the 
MW pole. They therefore appear close to the centre of the plot. In 
the bottom panel the most-favourable orientation of the spin of the 
DM halo of the MW is assumed. Here it is aligned with the average 
orbital pole of the MW satellite orbital poles from M08. Therefore, 
the histogram shows the distance of the MW satellite orbital poles 
measured from Galactic coordinates {l,b) = (f77.0° , — 9.4°). The 
MW satellite orbital poles cluster much more than expected from 
the Aquarius simulations. Half of the poles are found in the bin 
closest to the assumed halo spin, whereas the simulations expect 
only about fO per cent. All probability density distributions are 
normalised such that the sum over the area of all angular bins 
yields a probability of f . 



also show an over-abundance of counter-orbiting subhalos 
close to COS0H S = —1. In addition to the original data, 
in Figure [5] we have added the orbital poles of the MW 
satellites. The upper panel assumes that the spin vector of 
the DM halo of the MW is parallel to the MW disc spin, 
pointing to the galactic south pole. The bottom panel as- 
sumes the halo spin coincides with the average orbital pole 
of the six best-aligned MW satellites. It constitutes the most 
favourable orientation of the halo to explain the preferred 
orbital direction of the MW satellites from the simulations. 

In both cases, the observed distribution differs com- 
pletely from the models. In the first case, the MW satellite 
orbital poles peak close to cosOh s = 0, the MW equator. 
This is where the model distribution approaches its mini- 
mum, in contradiction to the observed distribution. In the 
second case, the observed orbital poles peak at the bin clos- 
est to the main halo spin by construction, which is also the 
maximum in the model distributions. But the peak contains 
50 per cent of the orbital poles and is therefore much higher 
than all simulated models, which predict only about 10 per 
cent in this bin. 

This illustrates that the observed situation does not fol- 
low the model. However, it does not rule out the possibility 
that the observed orbital poles are drawn from one of the 
model distributions. The question we need a quantitative 
answer for is: How probable is this arrangement? 



2.3 Models 

This section describes the different models whose angular 
momentum direction distributions will be analysed using the 
method presented in the next Section. For that, it is neces- 
sary to repeatedly draw eight angular momentum directions 
randomly from their distributions. 



2.3.1 Isotropic distribution 

Assuming an isotropic distribution of angular momentum di- 
rections for satellite galaxies, 8 random directions are drawn 
from an uniform distribution on the sphere. These are used 
as angular momentum vectors. 



2.3.2 Aquarius Simulations 

To test the results of Lll derived from the Aquarius cosmo- 
logical simulations, the angular momentum directions are 
drawn from the distribution of subhalo angular moment4f|. 
These are shown in their figure 4. Drawing from the whole 
population of subhalos and not only from the subhalos with 



^ The angular momentum data for model Aq-B2 were provided 
in a different orientation to those shown in figure 4 in Lll, there- 
fore the precise position of the main halo spin is not known. It 
was estimated by rotating the provided data so that it visually 
resembles the orientation in figure 4 of Lll. The distribution of 
angular momenta of this model in the re-production of their figure 
1 also closely follows the original of figure 1 of Lll, emphasising 
that the alignment is close to the correct one. Nevertheless, the 
results from this model are slightly more uncertain than those of 
the others. Fortunately, model Aq-B2 is the most isotropic one, 
so that this uncertainty has a negligible impact on the results. 
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the most massive progenitors should not pose a problem for 
this analysis, as Lll state that the latter trace the same 
structure of the whole subhalo sample. 

2.3.3 Via Lactea Simulations 

To provide a more complete picture of cosmological models, 
the subhalo angular momen tum directions derived from the 
two Via Lactea simulations dPiemand et al ] |2007l . l2008l . D07 
and DOS hereafter) are analysed, too. These are cosmologi- 
cal n-body simulations of MW-sized DM halos with a quiet 
merger history. They are therefore possible hosts of a MW- 
type galaxy in a ACDM-universe. The analysis is based on 
the data freely available on the website of the Via Lactea 
projeclQ for the two models Via Lactea 1 (VL-1, D07) and 
Via Lactea 2 (VL-2, DOS). The subhalo angular momentum 
directions needed for the analysis are then determined for all 
subhalos in the available data set that are within the virial 
radius. 

The halo spin is determined from the provided ran- 
dom sample of 10^ particles at redshift zero. The sum of 
the angular momentum of all particles within 0.25 ryjr (VL- 
1: Ty ir = 3S9 kpc; VL-2: r^ir = 402 kpc, see iKuhlen et al.l 
|2G08| ) is calculated, assuming each particle to have a mass 
such that the total mass of the simulation is reproduced. 
Then the sum of the subhalo angular momenta within the 
same radius is subtracted from the average particle angular 
momentum, which only leads to a small correction in the av- 
erage spin direction of l.S (VL-1) and 11.4 (VL-2) degrees. 
The direction of the resulting vector is then adopted as the 
(inner) halo spin direction, which is more closely aligned 
with the disc galaxy spin as discussed in Sect. 12. Il This pro- 
cedure is less sophisticated than the analysis Lll use, but 
can be easily done with the freely available data of the Via 
Lactea simulations and should give an estimate of the spin 
direction of the main DM halo. 

2.3.4 Tidal Models 

In addition to the cosmological simulations, models of a tidal 
origin for the MW satellites presented by Pll are anal- 
ysed. In these models, a target disc galaxy collides with 
a perpendicularly oriented infalling disc galaxy on a polar 
orbit. Pll included equal-mass and 4-to-l mass ratios for 
target-to-infalling galaxy. During the interaction, material 
is stripped off from the infalling galaxy and forms a disc 
of debris around the target galaxy, within which TDGs can 
form. Two populations of tidal debris form in most models, 
having pro- and retrograde orbits with respect to the orbit 
of the infalling galaxy. The models therefore not only pro- 
duce orbits similar to the majority of those described by the 
satellite galaxy orbital poles, but can also account for the 
counter-orbiting direction of the Sculptor dwarf galaxy. 

The distributions of angular momentum directions for 
these models are drawn from the tidal material, so from 
the angular momentum directions of individual particles. To 
show the evolution of the angular momentum distributions, 
three time-steps representing 5.0, 7.5 and 10.0 Gyr after the 
beginning of the model calculations are considered. These 

^ http://www.ucoIick.org/ diemand/vl/data.html 



are chosen to be well after the initial perigalactic passage of 
the two galaxies and, in the merger case, also after the final 
collision that happens after 2.5 to 4 Gyr. To include their 
angular momentum direction in the analysis, the particles 
are required to have a distance of at least 30 kpc from the 
central galajcy in order to avoid material from the galactic 
disc. If the model leads to a merger, particles with distances 
of up to 400 kpc are considered, while in a fly- by encounter 
particles can have distances from the target galaxy of up to 
half the distance between the two interacting galaxies. 

The results of all models from Pll are compiled in the 
appendix. In the following, four models will be discussed 
in detail. The four particular models have not been chosen 
to give the best results, but to show the typical range of 
fractions determined for the tidal models of Pll. They in- 
clude two 1:1 mass ratio models with the infalling galaxy 
oriented prograde, one which results in a fly- by ("l:l-fly by- 
pro" , called 5deg200vel in Pll) and one resulting in a merger 
("l:l-merger-pro", 7.5degl00vel in Pll) of the two colliding 
galaxies. In addition, two 4:1 mass ratio merger models are 
included. One of these ("4:l-merger-pro" , 7.5degl00vel in 
Pll), having a prograde infalling galaxy, is one of the mod- 
els with a good agreement. The other model ("4:l-merger- 
retro", lOdeglOOvel in Pll) has an infalling galaxy in ret- 
rograde orientation and results in the worst agreement with 
the MW orbital poles of all tidal models. The reason is that, 
while the model forms a very long tidal tail of more than 
300 kpc, it is dominated by a spheroid-like distribution of 
tidal particles out to about 50 kpc. 

2.4 Method 

To assess how likely it is to find a similar distribution of an- 
gular momentum directions as inferred from the eight MW 
satellite galaxies, it has to be estimated how often similar 
parameters can be produced by a given model. For this, a 
sub-sample of 8 different angular momentum directions is 
drawn from the distribution of angular momentum direc- 
tions in the model tested. 

Of the 8 known orbital poles of the MW satellites, 6 
are aligned with each other and the DoS normal direction. 
Analogously, for each possible combination of 6 out of the S 
angular momentum directions drawn from the models, the 
mean direction of the angular momenta is determined. Cen- 
tred on the se average directi ons, the spherical standard dis- 
tance Asph (jMetz et al.ll2007l ) of the 6 corresponding angular 
momentum directions is calculated. It is defined as 

/ELi [arccos ((n) ■ n,)]^ 
Aeph - Y , 

where hi are the angular momentum direction unit vectors, 
(h) is the unit vector pointing into their mean direction 
and '•' denotes the scalar produc t of the vectors. Note that 
the formula in lMetz et al.l (|2007l ) deals with axial data and 
therefore has to take the absolute value of the scalar product, 
in contrast to our case. Asph is a measure for the clustering 
of orbital poles around their average direction. The sample 
of 6 angular momentum directions leading to the smallest 
Asph is chosen. These are called the six best-aligned angular 
momenta from the sample of eight. The value for the average 
angular momentum direction and the spherical standard dis- 
tance of these best-aligned angular momenta is stored. The 
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spherical standard distance will be used for the clustering 
criterion in Sect. 13. Il 

Finally, the angular distance d of the average angular 
momentum direction of the six best-aligned angular mo- 
menta to the equator of the model is computed. In the 
isotropic case, the 'equator' is an arbitrary great circle as 
there is no preferred direction available. The great circle 90° 
away from the main halo spin is referred to as the 'equator' 
in the Aquarius and Via Lactea simulations, because the 
main halo spin aligns with the galactic disc spin, as was 
discussed in Sect. 12.11 In the case of the galaxy-interaction 
models, the 'equator' is the plane of the target galaxy which 
resembles the orientation of the MW disc. The angular dis- 
tance to the equator is used in the orientation criterion in 
Sect. 13.21 assessing whether the average angular momentum 
hints at polar orbits or not. 

When this algorithm is applied to the eight orbital 
poles of satellite galaxies of the MW, the same six satellite 
galaxies as reported by M08 are found to have the best- 
aligned orbital poles. Consequently, the algorithm gives the 
same parameters as reported by M08: A^]^ = 35.4° and 
dMW = 9.4°. 

This process of drawing 8 angular momentum directions 
and determining the parameters Agph and d for the best- 
aligned subsample of six of these is called one realisation. To 
determine the statistical properties, 10^ different realisations 
are produced for each model. The resulting distributions of 
the two parameters are shown in Figures [3] and [l] 

The method can be easily adjusted once more than eight 
satellite galaxy orbital poles become available through ob- 
servations of proper motions. The number of angular mo- 
menta drawn for one realisation (currently eight) would need 
to be increased, the number of best-aligned poles (currently 
six) might be changed, and the parameters A^JJ' — 35.4° 
and dMW = 9.4° adjusted to the observed situation. Note 
also that the analysis does not account for the alignment of 
the orbital poles with the DoS normal, nor does it factor 
in that the Sculptor satellite galaxy is counter-orbiting with 
respect to the six best-aligned orbital poles, but also orbit- 
ing in the DoS. If more satellite galaxies on counter-orbits 
are found, the analysis might need to be adjusted to account 
for this, e.g. by analysing not the directions of angular mo- 
menta, but the orbital axes (given by the angular momentum 
direction and its counter-direction), thus in effect combining 
CO- and counter-orbiting poles on one half-sphere. 



3 RESULTS 

The results are listed in Table [1] together with the total 
numbers of individual angular momentum directions of each 
model. In addition and for completeness, the results for all 
tidal models of Pll are compiled in the Appendix. 

3.1 Fulfilling the clustering criterion 

Figure [3] plots the resulting distributions of the spheri- 
cal standard distances Asph for the analysed models. The 
isotropic distribution peaks at about 55° and is spread be- 
tween 30° and 75°. 

The six Aquarius and the two Via Lactea simulations 
lead to Agph-distributions that are nearly indistinguishable 



from the isotropic distribution. The mean value of the spher- 
ical standard distances, Agph, are slightly lower, by up to 
about 3° compared to the isotropic case having A^ph = 54° . 

The first three tidal models, in contrast, show a com- 
pletely different behaviour in their Agph-distributions. They 
display a very strong peak at 5 to 15° in all three time steps. 
This peak gets lower for later time steps in all three models, 
showing that the angular momentum distributions widen 
with time. Their mean values Asph stay relatively constant 
at 31° (l:l-flyby-pro), 22° (l:l-merger-pro) and 18° - 25° 
(4:l-merger-pro). This shows that the tidal models result in 
much more concentrated distributions of orbital poles than 
the Aquarius and Via Lactea simulations. They are qualita- 
tively very different from the isotropic case and the concen- 
tration is stable for at least half a Hubble time. The Agph- 
distributions show tails spreading out to almost 90°. These 
arise from the fact that the angular momenta of the tidal 
debris cluster in two opposite directions, both being close 
to the central galaxy's equator. In addition to the strong 
clustering of orbital angular momenta, this pro- and retro- 
grade tidal debris can therefore explain the counter-orbiting 
direction of the Sculptor dwarf galaxy. However, this addi- 
tional aspect of the tidal models being consistent with the 
observed situation is not considered by the present analysis. 

The Asph-distribution of the last tidal model (4:1- 
merger-retro) approaches an isotropic-like shape for the last 
time steps. This illustrates that the spheroid-like distribu- 
tion of the merger remnant is dominating the angular mo- 
mentum directions. On average, the distribution is only min- 
imally more concentrated than those of the cosmological 
models. 

The orbital directions of satellite galaxies can be 
changed by a number of processes, like scattering with other 
satellite galaxies (whi ch m i ght e xplain the energetic or- 
bit of Sagittarius, see IZhadll998l ). precession due to non- 
axisymmetric potentials or tidal torques from neighbouring 
galaxies. For a clustered distribution of orbital poles, these 
processes lead to an increase of the spherical standard dis- 
tances. Therefore, the observed Agph-value of the MW satel- 
lites can be interpreted as an upper limit. It might have been 
smaller in the past, but not larger. 

To fulfil the clustering criterion, a realisation has to be 
at least as well clustered as the MW satellite orbital poles. 
The Asph-value of the MW satellites, 35.4°, is illustrated 
in Figure [3] by the vertical dashed line. All realisations with 
Asph ^ 35.4° are counted as passing the clustering criterion. 
Their number divided by the total number of realisations 
per model is /a, the fraction of realisations that fulfil the 
clustering criterion. The values of /a are compiled in Table[T] 
for all models shown in the plot. 

In the isotropic case, only 3 per cent of the realisations 
fulfil the clustering criterion. As can be expected from the 
longer tail towards lower Asph-values in the Aquarius sim- 
ulations, the fractions of realisations passing the clustering 
criterion are higher. However, this is not a strong effect, the 
increase compared to the value of isotropy is only a factor 
of 1.3 to 3.2. It might thus be concluded that filamentary 
accretion can account for a minor increase in the concentra- 
tion of angular momenta distributions of subhalos. The Via 
Lactea simulations lead to clustering results similar to those 
of the Aquarius simulations. 

Strongly clustered distributions of angular momentum 
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Table 1. Models and resulting fractions of realisations fulfilling the criteria 
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Mod el: Name of the mod el or s imulation. Reference : The original publication p resenting the respective model (Lll: UTovell et aLll201ll : 
D07: iDiemand et al.ll2007l : DOS: iDiemand et al.ll2008l : Pll: iPawlowski et al1l20l j) . N: Number of angular momentum directions to draw 
realisations from. Agph: Mean of the spherical standard distance distribution. Aspjj: Median of the spherical standard distance distribution. 
//^: Fraction of realisations fulfilling the clustering criterion (having a spherical standard distance of no more then 35.4°, the value for 
the MW satellites), f^: Fraction of realisations fulfilling the direction criterion (having an average angular momentum direction pointing 
no further away from the equator than 9.4°, the value for the MW satellites), /both- Fraction of realisations fulfilling both criteria, 
/indep = /a ■ fd- Fraction fulfilling both criteria if they were independent. All uncertainties are estimated assuming Poisson statistics. 

3.2 Fulfilling the orientation criterion 

The distribution of angular distances d of the average angu- 
lar momentum from the respective equator of the model are 
plotted in Figure|H The plot in fact shows sin(d), as this re- 
sults in bins representing rings of equal area measured from 
the equator. This gives a flat distribution for the isotropic 
case and therefore eases the comparison of the various dis- 
tributions. 

All Aquarius simulations have their maximum proba- 
bility density at the poles, far away from the equator. For 
the bins close to the equator their probability densities lie 
below the line corresponding to the isotropic case. They are 
therefore less likely to produce an average angular momen- 
tum vector close to the equator than the isotropic case. This 
results from the preferred alignment of subhalo orbits with 
the axis of the main halo spin, which has been illustrated in 
Figure [T] For model Aq-B2, its sin(d) distribution is almost 
flat, mostly following the isotropic case, but with a minor 
increase towards the poles. 

The Via Lactea models again have a behaviour similar 
to that of the Aquarius models, however the VL-2 average 
angular momentum vectors do not peak close to the pole, 
but shows a small increase towards the equator compared 
to the isotropic case. 

The tidal models again behave differently. Initially (5 
Gyr), the sin(d)-distributions of all four have their maxi- 
mum close to the equator (sin(ci) — 0). They widen slightly 
for the models l:l-flyby-pro and 4:l-merger-pro, the peaks 
become lower but do not move away from the equator much. 



directions arise naturally in tidal interactions. The first three 
tidal models easily fulfil the clustering criterion, for each 
time step more than half of all realisations have A^ph ^ 
35.4°: about 60 per cent in model l:l-fiyby-pro, 75 per cent 
in l;l-merger-pro and 80 per cent in the 4:l-merger-pro. 
There is only a small variation with time: /a tends to be- 
come lower with time (also compare the averages for all Pll- 
models in tables lAll and |A2)) . 



The last tidal model (4:l-merger-retro) starts with /a of 
34 per cent at 5 Gyr, but this value drops to only 8 per cent 
for the last time step, where it is comparable to /a of the 
cosmological models. This shows that the initially clustered 
distribution of angular momentum directions in this model 
quickly disperses because the spheroid-like component of the 
merger remnant dominates over the tidal tail material. If 
only the tidal tail particles are included in the analysis by 
demanding a minimum distance from the galactic centre of 
60 kpc, /a becomes 28 per cent for the final time step. This 
fraction is still relatively small compared to that of the other 
tidal models. 



The clustering criterion is completely independent of 
the orientation of the main galaxy. It alone already shows 
that reproducing the orbital-pole-distribution of the MW 
satellites in the Aquarius and Via Lactea simulations is very 
unlikely. 
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Figure 3. Distributions of the probability density (PD) of the spherical standard distance A^pij of the 6 most closely aligned angular 
momentum vectors out of 8 different, randomly drawn angular momentum vectors. In the first eight panels, the angular momentum 
vectors are drawn from the subhalo angular momentum vectors of the Aquarius simulations presented by Lll and from the Via Lactea 
simulations by D07 and DOS. The last four panels show the distribution that results when drawing angular momentum vectors from the 
particles in calculations of galaxy interactions by Pll. They each contain three distributions showing the time-evolution of each model. 
The thin, light gray histogram illustrates the situation at 5 Gyr after the start of the calculation, the dark gray histogram at 7.5 Gyr 
and the thick black histogram at 10 Gyr. The case of an isotropic sample of angular momenta is included as the thin black line in all 
panels for comparison. The spherical standard distance for the six best-fitting orbital poles in the MW, 35.4°, is indicated by the vertical 
dashed lines. To fulfil the clustering criterion, a realisation has to fall to the left of this line, this allowed region being highlighted by 
a shading. The distributions resulting from the cosmological simulations are nearly indistinguishable from the isotropic case. The vast 
majority of determined spherical standard distances is significantly larger than the observed one. This is completely different for the most 
examples of tidal material, which show very strongly concentrated distributions of angular momentum directions, in agreement with the 
distribution inferred from the MW satellite galaxies. 



This is different for model l:l-merger-pro, in which the dis- 
tribution peaks at sin(d) ~ 0.13 initially (5 Gyr), but moves 
to sin(d) « 0.25 (d = 14.5°) thereafter (7.5 and 10 Gyr), 
indicating the precession of the angular momenta of the par- 
ticles. The distributions of these three models drop to zero 
for sin(d) > 0.5 (d > 30°) at all time steps. The directions of 
their average angular momentum vectors thus all align close 
to the equator of the models. This is different for model 
4:l-merger-retro, which has sin(d)-values spreading over the 
whole range. Its maximum close to the equator is the lowest 
initially and the sin(d)-distribution of this model approaches 
the isotropic case for the last time step. 

The orientation of the average angular momentum vec- 
tor of the 6 best-aligned angular momentum vectors should 
be close to the equator in order to resemble the case of the 
polar orbits of the MW satellites, which have an angular 



distance of the average orbital pole from the MW equator 
of 9.4°, represented as dashed vertical lines in Figure [4] To 
fulfil the orientation criterion, the average angular distance 
from the equator has to be lower than or equal to this value. 
The fraction fd fulfilling this criterion is calculated similar 
to /a. 

The fractions of particles fulfilling the orientation cri- 
terion are compiled in Table [T] In the isotropic case, 16 per 
cent of the realisations are close enough to the equator. As 
expected, all Aquarius simulations and Via Lactea 1 show 
lower fractions, by a factor of 1.1 (Aq-B2) up to 2.7 (Aq-A2) 
less. Only Via Lactea 2 gives a slightly higher fraction, by a 
factor of 1.2. 

The tidal models again give the highest fractions of re- 
alisations fulfilling the orientation criterion. In model 1:1- 
fiyby-pro almost all (84-91 per cent) of the realisations lead 
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Figure 4. Distributions of the probability density (PD) of the angular distances of the average angular momentum vector from the 
equator for the 6 most closely angular momentum vectors out of 8 different, randomly drawn angular momentum vectors. Same panels 
as in Figure [S] The vertical dashed line indicates the distance of the average orbital pole of the MW satellites from the equator of the 
MW, d = 9.4°. To fulfill the orientation criterion, a realisation has to fall to the left of this line. As in Figure [T] the Aquarius and Via 
Lactea simulations produce results that make it unlikely to find a situation resembling the MW satellite orbital poles. Their subhalo 
angular momenta preferentially align with the main halos spin, which is perpendicular to the equator. An exception is the Via Lactea 2 
simulations, which shows a slight excess at low d. All other cosmological models have a lower number of average angular momenta close 
to the equator than the isotropic case. Most models of galaxy interactions however naturally come up with average angular momenta of 
their tidal debris close to the equator, resembling the observed case of the MW satellites. 



to average angular momentum directions close to the equa- 
tor. In model l:l-merger-pro this is only 15-60 per cent, 
because the peak of its d distribution moves further away 
from the equator for later time steps. For model 4:l-merger- 
pro the fraction is again high (53-83 per cent). In general, 
the fraction drops for later time steps (see also Tables lAll 
and lA2|l . The tidal model 4:l-merger-retro starts with a rel- 
atively high fraction at 5 Gyr (36 per cent), but this value 
drops quickly, such that at the final time step it is only 13 
per cent, comparable to the cosmological models and lower 
than the isotropic case. Considering only particles having a 
distance of at least 60 kpc from the galactic centre gives a 
higher fraction of 29 per cent at the final time step. This is 
again at the lower end of the range of values for tidal models. 

Doubling the accepted distance from the equator from 
rfMW ~ 9.4° to d = 18.8° increases fd in the merger case, 
leading to 75-92 (l:l-merger-pro) and 81-96 per cent (4:1- 
merger-retro). It about doubles fd for the the Aquarius and 
Via Lactes simulations, the isotropic case and the 4:l-mass 



ratio retrograde merger model, and increases by up to 10 
percentage points in the fiy-by case. 

3.3 Combining the criteria 

To resemble the distribution of MW satellite orbital poles, 
both criteria presented before have to be met simultaneously. 
If the direction of the average angular momentum and the 
spherical standard distance would be independent, the com- 
bined fraction should simply be the product of the two frac- 
tions. If the underlying distribution of angular momentum 
directions is non-isotropic, this will not be the case any more. 
Both the actual fraction of realisations that fulfil both crite- 
ria, fboth, as well as the expected fraction if the two criteria 
were completely independent, ftndep, are compiled in Table 

m 

As can be expected, in the isotropic model /both and 
/indcp are consistent with each other, 0.5 per cent of the re- 
alisations give parameters that fulfil both criteria together. 
It is highly unlikely to find the observed distribution of the 
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MW satellite orbital poles if the satellite orbits were dis- 
tributed isotropically around the MW. 

It has been claimed by Lll that cosmological simula- 
tions naturally give rise to "[...] distributions of coherently 
rotating satellites, such as those inferred in the Milky Way 
[...]". We tested this claim by determining /both, finding it 
to be falsified. All six Aquarius models result in /both-values 
that are below the one for the isotropic case. The model dis- 
playing the strongest alignment of subhalo angular momenta 
with the main halo spin, Aq-A2, results in a MW-satellite- 
like distribution of orbital poles in only 0.1 per cent of the 
realisations. Averaging over the six Aquarius models gives 
a likelihood of 0.33 per cent. VL-1 is consistent with the 
Aquarius models, having /both = 0.4. Only VL-2 is more 
likely to reproduce the observed situation than the isotropic 
case, but still only in 1.5 per cent of the cases. The reason 
is that our crudely derived halo spin direction (the sum of 
angular momenta of a subset of particles from the simula- 
tion) does not point into the direction in which most subhalo 
angular momenta point. 

Most tidal models, in contrast, result in much higher 
fractions passing both criteria. In the equal-mass fiy-by 
model (l:l-fiyby-pro), /both is as high as 60 per cent for the 
first time step. This value is only slightly reduced for the 
last time step, where it is 57 per cent. Therefore, this model 
will produce, insofar as we tested, angular momenta distri- 
butions with properties similar to that of the MW satellite 
orbital poles in most of the cases. The equal-mass merger 
case gives /both = 48 per cent for the first time step (5 Gyr) 
and 12 per cent for the last (10 Gyr). While lower than 
in the fiy-by case, this fraction is not low enough to rule 
out the model and is a factor of 27 larger than that of the 
best Aquarius simulation in this analysis and 8 times larger 
than the VL-2 result. The third tidal model (4:l-merger-pro) 
again gives very high /both-values, also dropping with time 
from 75 per cent at 5 Gyr to 46 per cent at 10 Gyr. 

Only the last tidal model (4:l-merger-retro) shows low 
/both-values, dropping from 19 per cent (5 Gyr) to only 1 
per cent (10 Gyr), a value still larger than (but only by a 
factor 2 to 3) the average of the cosmological models. If, 
however, angular momenta of the particles belonging to the 
spheroid-like component of the merger remnant are removed 
by demanding a minimum distance of 60 kpc from the cen- 
tral galaxy, then this model still has a /both of 13 per cent 
at the final time step. 

As discussed in Sect. 13.21 doubling the allowed distance 
from the equator increases fd and consequently also in- 
creases /both- This leads to /both = 67 per cent at the fi- 
nal time step for model l:l-merger-pro and to /both = 68 
for model 4:l-merger-pro. It only slightly changes /both for 
the other models, increasing the fraction by one percentage 
point in the fiy-by model and doubling the fractions for the 
Aquarius and Via Lactea simulations, the isotropic case and 
the last tidal model (4:l-merger-retro). 

Tables lAll and IA2I in the Appendix compile the results 
for all tidal models of Pll. High /both are common, some 
models even show values as high as 90 per cent or more. The 
/both-averages for the different model-types slightly decrease 
with time, but the fractions can increase for some models. 
Overall, they remain consistent for the different time steps. 
The /both-averages are highest in fly-by encounters (67 to 
52 per cent for the investigated equal-mass and 71 to 41 per 



cent for the 4:1 mass-ratio models at time steps 5 and 10 
Gyr, respectively). Retrograde mergers, on average, produce 
lower /both (31 to 33 per cent and 20 to 11 per cent) than 
prograde mergers (41 to 36 per cent and 44 to 27 per cent 
). The statistical properties of the distributions of angular 
momentum directions in the tidal models are stable over 
many Gyr. In order to reproduce the observed orbital-pole 
distribution, it is not necessary to fine-tune the point in time 
after the galaxy-galaxy interaction. 

Taken together, our investigations show that tidal mod- 
els are capable of naturally accounting for the observed dis- 
tribution of MW satellite orbital poles. In contrast to that, 
it can most certainly not be claimed that the Aquarius or 
Via Lactea models naturally produce the phase-space distri- 
bution of the MW satellites. There is no evidence in those 
ACDM models that filamentary accretion can account for 
the peculiar properties of the MW satellite orbits. Even an 
isotropic distribution of satellite angular momentum direc- 
tions has a higher likelihood of accounting for the MW dSph 
phase-space correlation than all of the Aquarius simulations. 
Tidal models, however, can reproduce the observed proper- 
ties with ease. 

Are there structures in the Aquarius simulations that 
resemble the MW satellite system? An answer might be 
sought in the spacial distribution of the DM subhalos, as 
shown in figure 6 of Lll. This figure depicts the positions of 
all subhalos, but only those are coloured which have orbital 
angular momenta aligned close to the axis of the main-halo 
spin (|cos(0_r/.s)| > 0.9). As these subhalos all move within 
approximately the same plane, they form a 'quasi-planar' 
distribution. However, it must not be forgotten that halos 
on these orbits are only slightly over-abundant, at most by a 
factor of about 2.3 compared to an isotropic distribution of 
subhalo angular momentum directions (for Aq-A2, see Fig- 
ure [TJ. The remaining subhalos are the majority and do not 
orbit in this plane. Thus, ignoring the suggestive colouring 
in figure 6 of Lll, it is obvious that the spacial distribution 
of subhalos does not resemble that of a disc of satellites but 
is much wider. The 'quasi-planar' distribution of Lll is, in 
fact, only a subset of the rather 'quasi-spheroidal' distribu- 
tion of DM subhalos in the Aquarius simulations. There is 
no mechanism to let luminous satellite galaxies only form in 
subhalos selected by their final position around the host halo 
they end up in. The chance to reproduce the observed, well 
pronounced DoS within the Aquarius models must therefore 
be extremely small. In addition, as discussed in Sections l2.H - 
12.21 the coloured subhalos of Lll emphasise a plane most 
likely oriented perpendicular to the DoS. 



3.4 A different equator 

Relaxing the demand of the orientation criterion, one might 
consider the hypothetical situation that the disc galaxy axis 
in the DM halo does not align with the main halo spin. 
The most favourable situation for the cosmological simula- 
tions would then be the case in which the equator lies where 
most subhalo angular momentum directions cluster. There- 
fore, the equator can be defined as the great circle connect- 
ing the main halo spin direction and its anti-pole with the 
direction in which the highest density of subhalo angular 
momentum directions is found (see figure 4 in Lll). Would 
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the Aquarius models result in satisfactory /both-values in 
this case? 

Even in this unlikely and extremely contrived situation 
for the cosmological models, their /both values can not com- 
pete with those of the tidal models. The highest value for 
/both, 3.6 per cent, is found for Aq-C2. The other models 
produce lower fractions: 2.7 (Aq-A2), 2.0 (Aq-E2), 1.7 (Aq- 
D2), 1.5 (Aq-F2) and 0.9 (Aq-B2) per cent. Averaged over 
the six models, the fraction is increased by only a factor of 
six, from 0.33 to 2.1 per cent. 



4 CONCLUSIONS 

A new method has been presented which estimates the likeli- 
hood of a given distribution of angular momentum directions 
to reproduce the peculiar properties of the observed MW 
satellite orbital poles, namely their preferred co-rotation on 
a near-polar orbit. For sets drawn from modelled angular 
momentum direction distributions, it tests both the close- 
ness of the average direction of the six best-aligned angular 
momenta to the equator and their spherical standard dis- 
tance. The method is easily applied to angular momentum 
distributions derived from any models and can be adjusted 
as more MW satellite galaxy orbital poles become available 
through observations of proper motions. 

If the MW satellite orbital poles were drawn from the 
orbital angular momenta of ACDM subhalos determined 
from the Aquarius simulations, then the substantial simula- 
tion effort has not been able to arrive at a distribution with 
similar parameters as the one observed around the MW. 
It was found that for the most likely orientation of the disc 
galaxy spin, a similar alignment has a probability of at most 
0.44 per cent to occur in the Aquarius simulations. This re- 
sult does not agree with the suggestion put forward by Lll 
that coherently rotating quasi-planar distributions of satel- 
lites such as inferred observationally for the MW arise nat- 
urally in simulations of a ACDM universe. 

The Lll cosmological models preferentially lead to 
slight over-abundances of angular momenta of subhalos 
aligned with the axis of the main-halo spin and so most 
probably with the axis of rotation of the central galaxy. This 
alignment is in stark contrast to the polar orientation of the 
orbits of the MW satellites and the polar orientation of the 
Disc of Satellites. Even the highly unlikely and contrived 
case of a perpendicular orientation of main halo and MW 
disc spin increases the probability to find the observed or- 
bital pole distribution by only a factor of 6 to about 2 per 
cent. 

An analysis of the two Via Lactea simulations gives sim- 
ilar results. With 1.5 per cent the VL-2 model has the high- 
est probability, as the sub halo spin directions in this model 
do not preferentially align with our crudely derived direction 
of the main halo spin. 

This result has consequences. iKeller et all (|2012l ) anal- 
ysed the distribution of globular clusters of the MW. Those 
globular clusters which are classified as young halo objects 
and are thought to have been stripped from accreted dwarf 
galaxies populate the same plane as the disc of satellite 
galaxies. I Keller et al.l (|2012l '). referring to the Lll results, 
have interpreted this as a sign of filamentary origin. They 
state: "Simulations have shown that the planar arrange- 



ment of satellites can arise as filaments of the surrounding 
large scale structure feed into the Milky Way's potential. 
We therefore propose that our results are direct observa- 
tional evidence for the accreted origin of the outer young 
halo globular cluster population". In the light of our results, 
the opposite is the case. The addition of more objects to 
the DoS puts even stricter demands on the models trying 
to explain the observed situation. As filamentary accretion 
already fails at explaining the orbital poles of eight MW 
satellite galaxies, more objects distributed in the same pla- 
nar structure (and therefore also orbiting in it, otherwise the 
structure would not be stable) only makes things worse. 

In contrast to a cosmological origin of the MW satellites, 
a tidal origin can easily reproduce the observed parameters. 
This was tested with models of interacting galaxies from Pll 
having different mass ratios and resulting in both fly-by and 
merger cases. Most of these models not only naturally pro- 
duce a clustering of the angular momenta, they also lead to 
discs of tidal debris around the central galaxy, similar to the 
DoS. Thus, a tidal interaction likely is the origin of the MW 
disc of satellites. The orientation of the orbital and spatial 
distribution of the tidal debris depend on the interaction 
geometry, which can well result in polar orbits and distri- 
butions. This is in agreement with the vast polar structure 
(VPOS) around the MW, which contains stellar and gaseous 
streams, globular clusters and satellite galaxies and their or- 
bits (Pawlowski et al. 2012). 

Indeed, other extra-galactic post-merger cases of 
aligned dSph satellites with tidal arms have emerged 
(Malphrus et al.'l99j: lBournaud et~alll2007l : [Galianni et all 
2010; Due et al. 201 ll ). strengt hening the notion that dSp h 
and dE satellites are old TDGs (|Okazaki fc Taniguchill2000l '). 
The implications of this for c osmological th e ory a n d funda- 
nicntal physics is discussed in lKroupa et al.l l|2010h : [Kroupal 
(|2012D . 

The percentages calculated in this work negate the sug- 
gestion that a satellite galaxy distribution like that of the 
MW arises naturally in a ACDM universe. Lll argue that 
the six Aquarius halos can be considered to be approxi- 
mately representative of the population of MW-sized ha- 
los as a whole. The similar results for the independent Via 
Lactea models seem to support this claim. Therefore, it can 
be expected that further simulations in the ACDM frame- 
work will not lead to significantly different results. Since 
the structure formation and merger history is similar in a 
warm-dark-matter cosmology (only the number of subha- 
los drops), the sam e conclusion holds true there as well. 
iKnebe et al.l l|2008l ) have compared the anisotropic spatial 
distribution of DM subhalos within their host halos for cold 
dark matter (CDM) and warm dark matter (WDM) simu- 
lations. They state that the spatial anisotropy of subhalos 
in the WDM model is in fact less pronounced than in the 
CDM case. Thus, CDM as well as WDM models appear to 
be ruled out. 

The tidal scenario, in which the MW satellite galax- 
ies are TDGs formed in a galaxy-encounter, is a promising 
alternative for their origin and able to ex plain the phase- 
space-corre lation of the DoS and VPOS jPawlowski et al] 
I201ll . l2012l ). 
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APPENDIX A: TABLES FOR TIDAL MODELS 

Tables \KT\ and IX2l compile the results of analysing three dif- 
ferent time steps of all 72 tidal-interaction models of Pll 
with the method presented in this paper. See Pll for de- 
tailed descriptions of the model set-up. The 4:1 mass ratio 
models have a target galaxy more massive than the infalling 
galaxy. In all encounters, the infalling galaxy is oriented per- 
pendicular to the target galaxy and approaches it on a polar 
orbit. The model names consist of two parts. The first, (e.g. 
'5deg') describes the angle between the velocity vector and 
the position vector of the two galaxies in degrees (here 5°), 
larger numbers give larger perigalactica. The second part 
(e.g. 'lOOvel') describes the initial relative velocity of the 
two galaxies with respect to the parabolic velocity Wparab in 
per cent (here 1.0 Wparab)- The fractions compiled in the ta- 
bles have been calculated from lO** realisations. Assuming 
Poission statistics, the uncertainties A/ can be calculated 
from the fractions / as A/ = y/J/10, where both A/ and / 
are given in per cent. 
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Time step 5 Gyr Time step 7.5 Gyr Time step 10 Gyr 



Type 


Model 


N 


/a[%] 


fd[%] 


/both[%] 


N 


/a[%] 


fd[%] 


/both[%] 


N 


/a[%] 




/both[%] 


fb 


2deg200vel 


1045 


99.4 


99.9 


99.3 


1136 


98.4 


99.3 


98.0 


1084 


98.3 


98.1 


97.0 


fb 


4degl75vel 


11882 


66.0 


79.9 


59.8 


12389 


29.5 


63.6 


26.5 


11793 


31.8 


49.8 


24.7 


fb 


4deg200vel 


14127 


34.7 


77.6 


33.8 


13732 


36.3 


70.8 


33.2 


12205 


39.9 


69.2 


36.4 


fb 


4deg225vel 


16477 


28.8 


83.9 


28.6 


16812 


27.4 


71.0 


25.7 


16276 


25.6 


56.9 


21.0 


fb 


4deg250vel 


11858 


38.3 


87.0 


38.1 


12016 


48.9 


92.2 


48.9 


11330 


45.0 


85.3 


444 


fb 


6degl75vel 


6284 


97.6 


86.7 


85.6 


7730 


42.3 


67.5 


37.1 


8245 


27.4 


41.9 


18.3 


fb 


6deg200vel 


6205 


97.9 


95.2 


93.7 


7730 


38.8 


65.0 


33.0 


7337 


31.5 


56.4 


24.9 


fb 


6deg225vel 


8659 


98.7 


94.7 


93.8 


8213 


85.4 


82.1 


73.5 


8509 


57.4 


60.1 


40.7 


fb 


6deg250vel 


4970 


94.3 


96.0 


91.4 


5067 


77.4 


62.0 


50.3 


no data, 


galaxies left model volume 


fb 


8degl75vel 


4522 


96.2 


93.1 


91.0 


5129 


80.9 


84.3 


73.5 


4815 


74.2 


76.5 


63.4 


Average 




75.2 


89.4 


71.5 




56.5 


75.8 


50.0 




47.9 


66.0 


41.2 


mp 


2.5deg050vel 


66917 


32.8 


12.7 


6.9 


62262 


23.9 


27.4 


10.8 


63893 


14.8 


28.2 


6.2 


mp 


2.5degl00vel 


81684 


44.9 


23.9 


16.2 


82953 


31.2 


17.7 


7.4 


81788 


25.9 


27.2 


7.5 


mp 


5deg050vel 


64938 


38.4 


24.9 


16.6 


60115 


29.0 


17.8 


7.9 


64109 


24.3 


28.9 


6.9 


mp 


SdeglOOvel 


86824 


86.1 


68.5 


63.5 


82448 


77.8 


58.5 


51.1 


81435 


75.9 


59.9 


50.4 


mp 


7.5degl00vel 


93940 


82.6 


82.3 


74-6 


80454 


89.1 


74.6 


70.9 


77768 


74.0 


53.9 


47.0 


mp 


lOdeglOOvel 


83850 


89.8 


90.0 


85.7 


72029 


77.1 


74.3 


64.0 


68751 


65.2 


53.2 


43.4 


Average 




62.4 


50.4 


43.9 




54.7 


45.0 


35.3 




46.7 


41.9 


26.9 


mr 


2.5deg050vel 


69336 


27.8 


72.8 


25.9 


65802 


25.7 


45.9 


18.3 


72171 


25.7 


43.1 


17.2 


mr 


2.5degl00vel 


80994 


42.1 


53.0 


29.1 


78535 


21.6 


41.7 


13.4 


78556 


22.9 


45.2 


15.5 


mr 


5deg050vel 


64109 


30.8 


62.8 


25.3 


61710 


25.3 


46.7 


16.9 


64159 


20.9 


45.3 


14.5 


mr 


SdeglOOvel 


72679 


22.9 


42.8 


15.7 


61278 


19.7 


34.2 


10.6 


59795 


23.1 


41.1 


15.3 


mr 


7.5degl00vel 


70490 


56.5 


4.3 


3.3 


51645 


41.1 


36.3 


18.4 


49391 


17.4 


19.4 


4.4 


mr 


lOdeglOOvel 


75431 


32.8 


36.2 


18.8 


60248 


16.2 


28.4 


6.9 


61157 


7.9 


13.8 


0.9 


Average 




35.5 


45.3 


19.7 




24.9 


38.9 


U.l 




19.6 


34.6 


11.3 



Type: Type of interaction. Either a fly- by with the infalling disc in prograde orientation (fb), or a merger with the infalling disc in prograde 
(mp) or retrograde (mr) orientation. Model: Name of the model or simulation. To show the evolution of the models, the following data 
is given for three different time steps in each model, except when stated otherwise. A'^: Number of angular momentum directions to draw 
realisations from. //\: Fraction of realisations fulfilling the clustering criterion (having a spherical standard distance of no more then 35.4°, 
the value for the MW satellites) . fd ■ Fraction of realisations fulfilling the direction criterion (having an average angular momentum direction 
pointing no further away from the equator than 9.4°, the value for the MW satellites), fhoth- Fraction of realisations fulfilling both criteria. 
They are printed in italic to ease the comparison of the different columns. 
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Table A2. Equal-mass tidal models and resulting fractions of realisations fulfilling the criteria 









Time step 5 Gyr 






Time step 7.5 Gyr 






Time step 10 Gyr 




TVDG 


Model 


N 


/a[%] 


fd[%] 


J both L ' 


N 


/a [%] 


fd[%] 


fhr^tl, f%l 

J botn L ' 


N 


/a[%] 


fd[%] 


fhotl, f%l 

J botn L ' 


fb 


5degl80vel 


3279 


86.5 


96.6 


85.8 


3391 


87.0 


94.9 


85.4 


3355 


82.2 


89.2 


78.1 


fb 


5deg200vel 


5821 


59.4 


91.3 


59.4 


5841 


59.6 


86.6 


58.9 


5756 


58.7 


83.4 


57.5 


fb 


5deg220vel 


3991 


46.5 


90.4 


46.5 


3800 


55.0 


89.9 


54.9 


3642 


44.8 


82.0 


44-3 


fb 


5deg240vel 


2189 


76.5 


99.2 


76.5 


1871 


71.2 


97.3 


71.2 


1596 


61.9 


92.3 


61.8 


fb 


6degl80vel 


2108 


97.7 


99.1 


97.3 


1998 


93.7 


98.9 


93.4 


2125 


86.0 


95.6 


85.1 


fb 


6deg200vel 


4701 


76.3 


93.5 


75.6 


4671 


56.1 


81.1 


54.1 


4793 


44.1 


70.4 


41.3 


fb 


6deg220vel 


3732 


45.8 


91.8 


45.8 


3840 


40.9 


82.0 


40.6 


3805 


32.9 


71.8 


31.5 


fb 


6deg240vel 


3797 


31.0 


95.2 


31.0 


3464 


32.7 


92.8 


32.7 


3093 


37.2 


84.4 


36.7 


fb 


7degl80vel 


1636 


98.3 


99.7 


98.2 


1712 


93.6 


97.1 


92.5 


1959 


65.1 


70.5 


57.8 


fb 


7deg200vel 


3010 


68.0 


94.5 


67.8 


3361 


48.2 


77.7 


45.8 


3297 


42.3 


67.4 


37.8 


fb 


7deg220vel 


2209 


50.5 


95.6 


50.5 


2406 


35.2 


80.6 


34.7 


2277 


31.1 


71.8 


29.2 


fb 


7deg240vel 


2423 


39.8 


94.2 


39.8 


2212 


45.2 


94.6 


45.1 


2125 


57.5 


94.5 


57.3 


fb 


8degl80vel 


1353 


99.1 


99.2 


98.6 


1433 


95.5 


97.9 


94.6 


1666 


75.7 


82.3 


69.9 


fb 


8deg200vel 


1939 


82.5 


98.4 


82.4 


2077 


69.1 


88.8 


66.7 


2232 


43.6 


72.0 


40.6 


fb 


8deg220vel 


1257 


66.2 


96.2 


66.0 


1292 


63.8 


88.6 


62.5 


1214 


51.3 


75.6 


47.3 


fb 


8deg240vel 


1298 


42.3 


92.1 


42.3 


1190 


52.3 


96.3 


52.2 


1112 


61.9 


93.9 


60.9 


Average 




66.6 


95.4 


66.5 




62.4 


90.3 


61.6 




54.8 


81.1 


52.3 



mp 


0deg050vcl 


11212 


53.1 


2.3 


2.1 


16610 


42.3 


34.0 


12.1 


17331 


56.9 


53.0 


33.8 


mp 


OdeglOOvel 


14095 


51.3 


20.8 


18.0 


15821 


51.4 


11.9 


6.3 


16217 


43.6 


51.9 


27.9 


mp 


2.5deg050vel 


11452 


46.7 


1.4 


0.8 


17368 


39.3 


27.6 


7.6 


14000 


45.5 


37.9 


17.3 


mp 


2.5degl00vel 


11353 


30.3 


19.4 


10.0 


11652 


28.2 


16.9 


2.9 


12825 
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